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Finite temperature Green's function technique is used to calculate the energies and damping rates 
of elementary excitations of the homogeneous, dilute, spin-1 Bose gases below the Bose-Einstein 
condensation temperature both in the density and spin channels. For this purpose the self-consistent 
55 ' dynamical Hartree-Fock model is formulated, which takes into account the direct and exchange 

processes on equal footing by summing up certain classes of Feynman diagrams. The model is 
shown to fulfil the Goldstone theorem and to exhibit the hybridization of one-particle and collective 
excitations correctly. The results are applied to the gases of 23 Na and 87 Rb atoms. 

' PACS numbers: 03.75.Mn, 03.75.Kk, 03.75Hh 

o 

O \ I. INTRODUCTION 

The experimental realization of Bose-Einstein condensation with ultracold, dilute, alkali gases in magnetic traps 
[1, 2] made it possible, among other things, to meaningfully compare the results of experiment and theory. The reason 
lies in the diluteness of the gas sample, so the interaction between the atoms of the gas can be modeled by a week 
pseudopotential, in most of the cases, which can be considered as perturbation in calculations. With the development 
of optical traps it became possible to confine not just one or two of the magnetic states |/, m) of the atom (the 
f*^; ■ weak-field seeking states), but all of them for a given hyperfine spin / giving rise to a so called spinor condensate 
[3-7]. Because of these internal degrees of freedom such systems can have much wider dynamical properties than 
those confined in magnetic traps. The theoretical study of spinor Bose condensed systems was initiated by [8, 9]. 

There are two different phases of spin-1 systems in the presence of a Bose condensate distinguished according to the 
total magnetic moment of the system [8, 9]. Namely, the ferromagnetic phase, which is characterized by the existence 
of a macroscopic magnetic moment and the polar phase where the system has no macroscopic magnetic moment at 
all. In the absence of an external magnetic field the type of the phase of the spin-1 Bose system depends only on 
the interaction potential of the atoms, which can be accurately characterized by s-wave scattering lengths ao and ai 
responsible for scatterings in the total hyperfine spin channel zero and two, respectively [8]. If 02 > ao the system 
prefers the polar phase, while in the case of 02 < oq the system is ferromagnetic [8]. For both cases there exist species 
of alkali Bose atoms. For example the / = 1 hyperfine spin state of 23 Na realizes the polar phase, while that of 
87 Rb the ferromagnetic phase. Since in the absence of a magnetic field the phase of the system does not depend on 
experimentally tunable parameters, one can use the word 'case' instead of 'phase'. 

The elementary excitations of dilute Bose gases have been studied for a long time [10, 11] , since these gases represent 
the simplest prototypes of symmetry breaking quantum systems. The zero temperature elementary excitations of these 
systems are described reasonably well by the Bogoliubov Hamiltonian [10]. It is not straightforward to generalize 
Bogoliubov's theory to finite temperatures in a consistent way [12-19]. The simplest finite temperature approach is 
the self-consistent Hartree model, also known as the random phase approximation [13, 16, 19]. 

The inclusion of exchange processes is not trivial in a consistent way using the finite temperature Green's function 
technique as shown in [17] for a scalar Bose-gas. (See for a formulation of the Hartree-Fock theory for such systems 
giving the density response functions [20]) For instance self-consistency has proven to be essential for developing a 
theory satisfying general requirements. As a matter of fact applying unperturbed (i.e. free particle) propagators as 
internal line contributions leads to an instability for long wavelength fluctuations in case of repulsive interactions. 
(Note, however, that for attractive interactions such model has been used to describe the collapse of the gas phase 
[21].) 
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In the present paper we generalize the model presented in [17] with systems with scalar condensate to the spin-1 
Bose gas. Throughout the paper we will work in the orderd phase. We start from the self-consistent Hartree- 
Fock approximation for the equation of state, whose iteration we reveal which contributions are included in the 
self-consistent procedure. On this basis the self-energies of various one-particle Green's functions can be given in 
harmony with the equation of state. Furthermore the obtained self-energies fulfill the Goldstone theorem, which is 
the consequence of the breaking of the gauge and spin rotation symmetries. This ensures gapless excitation spectra. 
Having determined the one-particle propagators, the expressions of various correlation funcions will be given. The 
requirement of the coincidence of the one-particle and certain collective excitations is satisfied and the model resulted 
fits into the dielectric formalism. (For a detailed discussion of the dielectric formalism for spin-1 systems see [19] and 
for a general review [13].) One has to point out that the model developed can be regarded as an improved mean-field 
one and as such looses its validity near the phase transition point. 

Since the present work deals with homogeneous systems, the results are applicable to trapped gases only when a 
local density approximation is relevant or when the quantities in question are not sensitive to the finiteness of the 
system. In analogy with what proved to be the situation in the case of gases with a scalar condensate (see [17, 22] 
and references therein) one can hope that the damping of the elementary excitations belongs to the latter category. 
Concerning the local density approach it remains to be seen whether such experiment which provided the frequency 
of local excitations [23] can be carried out for spin-1 gases. Anyhow, we evaluate the characteristics of the excitations 
for parameter values which might be relevant to applications. Finally we note that a generalization of the model to 
finite systems is straightforward like for systems with a scalar condensate [22, 24, 25]. 

We hope to carry out the calculation directly to trapped systems in the future. At this point one has to emphasize 
again the self-consistent nature of the model which is very important for applications to confined gases. 

The outline of the paper is as follows. In Sec. II. the formulation of the Hartree-Fock approximation is given for 
the spin-1 Bose gas. It includes the specification of the Hamiltonian, the dealing with the symmetry breaking by a 
canonical transformation, the equation of state, the definitions and properties the Green's functions and of various 
correlation functions. In this section we specify what building blocks are to be used for the determination of the 
various quantities such that 

• being consistent with the implicit summation occurring in the Hartree-Fock propagator due to the self- 
consistency condition, 

• fullfilling the hybridization property. 

In Sec. III. the Hartree-Fock model is applied to a ferromagnetic system, the gas of 87 Rb atoms [Sec. Ill A], and to 
the 23 Na gas, which realizes the polar phase [Sec. IIIB]. Sec. IV. is devoted for further discussion also focusing on an 
overview of the physical properties of the various collective modes considered in Sec. III. Some details are relegated 
to the Appendix. 

II. FORMULATION 

We consider a system of spin-1 bosons in a unit volume box with periodic boundary conditions. For such a system 
the plausible basis set of one particle states to build the Fock-space from are plane wave states. The creation and 
destruction operators are a£ (k) and a r (k), which create and annihilate a one-particle state with momentum k and spin 
projection r, respectively. Throughout the paper Roman indices are used for distinguishing between spin components 
of a quantity. These indices will always take values from the set {+,0,—} referring to the eigenvalue of F z . Unless 
otherwise stated we understand automatic summation over repeated Roman indices. For our purposes, while the 
temperature of the system is low, and the gas is very dilute, the system can be well described by the following low 
energy, effective, grand-canonical Hamiltonian [8]: 

h = 5> k - M )4(kK(k) + i Yl 4'( k i)4(k 2 )v;';' s 'a s (k 3 ) as ,(k4), (i) 

k ki+k2=k 3 +k 4 

where ek = H 2 k 2 /2M is the kinetic energy of a particle with mass M, and fi is the chemical potential. We restrict 
ourselves to the Hartree-Fock approximation, and for such purposes the interatomic potential in Eq. (1) can be 
modeled by the momentum independent pseudopotential [12, 26]: 

V?/ = c n 5 rs 5 r . s > + c s (F) rs (F) rV , (2) 
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with parameters: 



Anh 2 a + 2a 2 
~M 3 : 
47rfi 2 a 2 — a 
~~M 3 ' 



(3a) 
(3b) 



The parameters a and a 2 are the scattering length in the total hyperfine spin channel zero and two, respectively. 
Note, that c s > for the polar case and c s < for the ferromagnetic one. In our representation the one-particle spin 
operators are 
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Let us denote the temperature with To, where the symmetric (uncondensed or high temperature) phase becomes 
unstable. To treat the system in the Bose condensed phase one can introduce the following new set of operators with 
a canonical transformation: 



6 r (k) = o r (k) - <5 k! ovA/VoCr, 
&t(k) = aj(k) - <5k,o75Vo£- 



(5a) 
(5b) 



The Nq and £ r parameters appearing in the transformation (5) describe the number and spinor of particles in the 
condensate, respectively. The spinor is normalized to unity. The spinor describing the polar state can be taken as 
C = (0, 1,0) T , while the ferromagnetic spinor as £ = (1,0, 0) T . Here and in the following, the 'T" in the superscript 
means transposition. The equation of state consists of two parts. The first one is the condition 



<Mk)) = (&t(k)) = o, 



(6a) 



which provides an equation among N and the chemical potential (see e.g. Ref. [19]). The other part is the expression 
of the particle number. In our case 



N 



(4(kK(k)) = N (T^) + ]T (bt(k)b r (k)) . 



(6b) 



where we have used Eqs. (5) to arrive to the right hand side of the above equation. 



A. Green's functions 



To calculate the grand-canonical averages in the equation of state (6) we use the method of finite temperature 
Green's function theory. The Green's functions of the system are defined by: 

^S(k,r):=-(T T [67(k,r)6f (k,0)]), (7) 

where r stands for the imaginary time and T T is the r ordering operator. The Greek indices are introduced for 
abbreviation, with &*(k) = 6 r (k) and b~ 1 (k) = &£(— k). Since they show the direction of propagation of the Green's 
function, we also will refer to them as direction indices. The convention of automatic summation is also understood 
for the Greek indices unless stated otherwise. 

These functions are periodic in r with period (5h. It is more convenient to use the Fourier series instead of the 
imaginary time dependent functions: 

^?(k,zu,„) = -Lj" K d T e^g r 4(KT), (8) 

with lo„ = 2mr/PH the Bose discrete Matsubara frequency [27]. 

The perturbation series of the Green's functions permit a simple partial summation. The Green's functions satisfy 
the generalized Dyson-Beliaev equations [13, 19, 27]: 



<?^(k, iu n ) = giZ hg (K iu n ) + S[ohp(k, i"n)s£' (k, iw„)0#(k, icj n ), 



(9) 
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where 

9 ^C.»J = CT^ < 10 > 

is the free propagator, and S^ CT S (k, ioj n ) is the self-energy, i.e. the sum of the contribution of the one particle 
irreducible Feynman diagrams, which connect to two external lines, one from the left, with indices r' , p, momentum 
k and frequency u) n , the other from the right, with indices s',a, momentum k and frequency u) n (we will call the 
one particle irreducible graphs simply as irreducible graphs). Note, that the chemical potential is included in the 
self energy instead of the frcc-proagator (see, e.g., Ref. [19]). Both the Green's functions and self-energies obey the 
following symmetry properties [19]: 

^S(k,iw„) - S£(k,tw n ) = g s _ r s ,_ 7 (K -iw„), (11) 

where E can also stand instead of Q in Eq. (11). Since spin and space variables are not coupled, the Green's functions 
and self-energies depend only on the modulus of the momentum. 

According to the usual procedure elementary excitations can be obtained. Continuing the Green's functions in their 
frequency variable from the imaginary Matsubara frequencies to the complex plane, requiring to be analytical on the 
upper-half plane, one arrives at the retarded Green's functions, the poles of which (in the lower half-plane) give the 
various one-particle elementary excitations of the system. 



B. The equation of state 



Now let us return to the problem of determining the equation of state (6) of the symmetry breaking dilute spin-1 
Bose gas in the Hartree-Fock approximation. Indeed, the equation of state in the Hartree-Fock approximation is the 
system of nonlinear equations resulting from the self-consistency condition of the Hartree-Fock propagators along with 
Eqs. (6). 




FIG. 1: The Feynman diagrams contributing to the self-energy S^f of the internal line (a) and the contribution of the tadpole 
diagrams (b). The heavy line depicts the internal line, while the thin line corresponds to the free propagator. The circle denotes 
a condensate operator, the dashed line is for the interaction and the triangle corresponds to —p. 



We denote the Hartree-Fock propagators (which will play the role of internal lines) by (/"(k, iu n ) and define it 
with the following equation: 



£?[i s (k,iu; n ) = 



iu n -h - En 



(12) 



where EJf is the momentum and frequency independent self-energy of the internal line. The Feynman diagrams 
contributing to the self-energy of the internal line can be seen in Fig. 1. The corresponding equation reads as: 



-ps rs + (N {r>{.> + H r ' s ') (v;/ + vz 



where H r s is the contribution of the loop (appearing in the Hartree and Fock terms): 



V f A 1 
lim / 

v-^oj (2tt) 3 fill 



£[-£!(q,«'n) 



(2^) 2 A 3 



r(3/2)F(3/2,/3fi£-), 



(13) 



(14) 
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with 

the thermal wavelength, T(s) the gamma function and F(s, 7) the Bose-Eintein integral with parameter s and argument 
7 [28]. 

The value of the chemical potential is determined by the requirement Eq. (6a), or equivalcntly by the vanishing of 
the contribution of tadpole diagrams [19]: 

= E' Q1 = ft" 1 V^bCr [ - rfrs + jVoCr'Cs'K"/ + H* (V/ + K''/) ] ■ (16) 

Since we work in the ordered phase, N ^ everywhere assumed. The simultaneous solution of equations (12), 
(13), (14), (16) (which is the consequence of Eq. (6a)), and Eq. (6b) provides the self-energy and propagator of the 
internal lines. 

For a ferromagnetic system the condensate spinor can be chosen to be ( r = (1,0, 0) T . According to Eq. (16) the 
chemical potential for this case reads as: 

H = c„ (N + 2H++ + H 00 + H-) + c s (N + 2H++ + H°°-H—), (17) 

while from Eqs. (13) and (17) the self-energies take the following form: 

h%ti = (c n + c s )No, (18a) 

fiE°° = -{c n + c s )H++ + (c - c s )H 00 + 2c,H—, (18b) 

fiSrf = ~(cn + 3c s )iJ++ + (c + 3c s )H— - 2c s N . (18c) 

The contributions of the single loops are also related to the self-energies through Eq. (14). There is a final equation 
from Eq. (6b) reading as: 

N = N + H++ + H 00 + H-. (19) 

The system of equations (17), (18), (14) and (19) form the equation of state in the following way. One can express the 
condensate density with the help of N, H ++ , H oa , H~~ from Eq. (19) and substitute it to Eqs. (18). These equations 
now form a closed set of equations toether with Eq. (14). One can solve them for H ++ , H 00 , H~~ as functions of 
T, N, and substitute them to Eq. (17) to arrive at the equation of state /U = /z(T, N). An alternate way is to express 
iVo from Eq. (17) and solve the equations for the loop terms with T, \x fixed, and then express N = N(T, fi). The two 
equations of state are necessarily equivalent. 

For a polar system the condensate spinor can be taken as ( r — (0,1, 0) T . The chemical potential for this phase, 
from Eq. (16), reads as: 

H = c„ (7V + H++ + 2H m + H-) + c s (H++ + H" ) . (20) 
The self-energies (13) in this case take the following form: 

= c n A^ , (21a) 
h%++ = (c - c s ){H++ - if 00 ) + c s N 0l (21b) 

where we have used the specific form of the chemical potential (20) and that E^ 4 " = EJV - and H ++ = H~~ . Equation 
(6b) can be expressed as: 

N = N + 2H ++ + H 00 . (22) 

Expressing No from Eq. (22) and substituting it to Eqs. (21) the latter ones can be solved together with Eq. (14) 
for H 00 and H ++ . They can be used in Eq. (20) to express the equation of state as \x = fj,(T, N). The alternate way 
is also possible here, i.e. to express the total number of particles as N = N(T,fj,). 



C. Effective interactions 



In this subsection we discuss the effective interactions which will play an important role in the following. The 
self consistent solution of the Hartree-Fock propagators contains a partial summation of infinite number of Feynman 
graphs. This partial summation is visible from the iterative solution of the Hartree-Fock propagators. With iteration 
one arrives to effective interactions which are the contributions of repeated exchange and direct terms. 
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1. Repeated exchange interactions 



Denote the contribution of multiple exchange interactions with the quantity 6£f s , (k\,k 2 ,q), defined by the Feynman- 
graphs in Fig. 2. For the sake of brevity one can introduce the four- vector notation: q = (q,iv n ) and k = (k, iu) n ). 



ki+q 



s I I r 

r' I s' 
«-l 1— « 



fcl+q fci 
+ 



fcl+q' fei+q 



&2 ^2+9 &2 ^2+9 &2 k 2 +q' k 2 +q 

FIG. 2: The Feynman diagrams defining ©"^(fci, fa, </)■ The shaded square depicts 0. 
The equation for reads as: 

Q r r f sl (k u k 2 ,q) = V r V s , + [ D q '\e a r f b (ki,k 27q ')p {0) b c a d (k 1+ q',k 2 +q')V^ 



where 



d 3 fc 

(2tt) 3 



and 



f(0)*^(Ai,fe) := -ft- 1 e?i ri (fci)e;r(fc 2 ) 



(23) 

(24a) 
(24b) 



were also introduced. The contribution of the bubble graph can be calculated from p( ), which we give here for future 
reference: 



/ 



0flP(o££(«,* + 9) = n (o ££(fc). 



(25) 



Due to the conservation of the spin z component Eq. (23) decouples to independent matrix equations. Let us define 
the following matrices for spin transfer: 



(26a) 
(26b) 
(26c) 



°§ = 


- e++ 

+± 


B 0+ 

^+0 




e# 


B 00 

e°~ 


©o-° 

e:: _ 


% = 


C n + 







Cs 




Cs 







c s 





J £(o) 







P( )oo c 
Pm 



For +1 spin transfer one defines 



u 00 ^0- 
+0 o00 



C s C n 



P(0)% 
P(0) - 



And finally, for +2 spin transfer one introduces 



Qy — r _ c 
Q £(o) =P(o)+^- 



(27a) 
(27b) 
(27c) 

(28a) 
(28b) 
(28c) 
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With the definitions (26), (27) and (28) one can cast Eq. (23) to the following integral equations for n G {0, +, Q}: 



Dq n Q.{k 1 ,k 2 ,q') 5(q-q')2_- ™P( ) ( fc i + l' , k 2 + q') ™Z 
The (5-function is normalized to unity with the definition 

It is straightforward to check that the solution of Eq. (29) reads as: 



l V. 



n @(k u k2,q)= n V- V \k2-h), 



with 



V(fc)= I- "S ( o)W n z 



(29) 



(30) 



(31a) 



(31b) 



The matrix n U_(o) is understood as the integral of the matrix ™P(o) according to Eq. (25). The matrices ™e* can 

be interpreted as dielectric functions with respect to repeated exchange interactions. It is important to note that 
™6(fci, q) does not depend on the amount of momentum transferred (q) at all, as is visible from Eq. (31a). 



2. Repeated direct interactions 

For the determination of the effective interaction representing multiple direct interactions, first consider the regular 
polarization functions, defined by the Feynman diagrams depicted in Fig. 3, or with the following equation (already 
separated by spin transfer): 



*JI W (fc)= "2( )(fc) + / DqDqi n g {n) {q,k + q) n g(q, k + q, q' - q) " g {0) {q' , k + q% (32) 



which leads to 

"n (r) (fc)= n i~\k) "n (0) (fc)- (33) 




■o- + --<n> 



FIG. 3: The Feynman diagrams defining the regular polarization H^J r . The hatched poligon represents the regular polarization 
function. 



The direct effective interaction is then defined by the Feynman diagrams in Fig. 4, or by the equation 

W r 7 s , (k) = V r 7 s , + W2 (k) n^" (k) V r % ■ (34) 

Equation (34) splits into separate matrix equations according to the amount of spin transferred. Introducing the 
quantities: 



°W 



wit w y- w++ 

Wl\ VF ° ° W°°_ 



Kt <- + 

wj w -_° 



Q C = 0, 



Cn H~ C s C n C n CS 
Cn C n C n 

c n c s c n c n ~\~ c s 



C.s C s 
C.<? Co 



(35a) 

(35b) 
(35c) 



FIG. 4: The Feynman diagrams defining the direct effective interaction. The wiggly line represents the direct effective interac- 
tion. 



the resulting equations are 

n K(k)= n g+ n K(k) n U (r) (k) n g (36) 

for n e {0, +, Q}. From Q C = it immediately follows that Q ]£ = 0. 
The solutions of Eq. (367~reads as: 

n K(k) = n g n g (rr \k), (37) 

with 

"i (r) = "I - ™fl (r) n Q ( 38 ) 

3. The equation of state in terms of the effective interactions 

The final structure of the Hartree-Fock self-energies and tadpole diagrams, as a result of the iteration of the Hartree- 
Fock propagator, are depicted in Fig. 5. a) and b), respectively, with the help of the effective interactions. For the sake 
of simplicity we have omitted the spin indices from the figure. Note, that in the final structure of the self-energies in 
Fig. 5. a) the Hartree-Fock propagators do not join vertices of the same interaction line, such propagators are always 
unperturbed ones. The iteration, of course, can be continued by substituting the remaining Hartree-Fock propagators 
with self-energies given in Fig. 1 a) in the self-energies, a procedure without end. These further iterations, however, 
do not generate new blocks, i.e. the building blocks remain the effective interactions and W . 

4- Anomalous vertex functions 

For the sake of a simplified notation it is convenient to introduce the anomalous vertex functions A^(fc), and 
A"",(fc) , whose general definitions can be found in [19]. In the Hartree-Fock approximation the Feynman diagrams 
for A^ r 1 are depicted in Fig. 6. The equations for the anomalous vertices in this approximation read as: 

Ki(k) = A£ )0il + J Dq p^ )cd (q, k + q)Gti(q, k + q, -q)A^' aA , with A£ )o>1 = 5 B , a C r y/No, (39a) 

K r ,-i(k) = A(o) a> _i + J Dq fft )cd (q, k + q)Q%{q, k + q,-k- q)A( Q ^_ 1 , with A^ )o ,_! = 5r,aC. V^b, (39b) 

K'l' (k) = A$ rV + J Dq A^„eS(0, k, q)p% r , s , (q, k + q), with A^ r , s , = S a , s> Cr' V^o, (39c) 

A^(k) = A^X'+ J D Q A (o)r\®T^ with A°$ B ,=8 a , r .C,'VNo- (39d) 

Note that in the diagrams for (39c) and (39d) the type of left and right connections are interchanged as compare to 
the diagrams for (39a) and (39b), respectively. 



Q Q 



-^*X*JL*x*j 



(a) 



+ 



+ 












+ 



QO 




+ 




(ML 



(b) 



FIG. 5: The Feynman diagrams contributing to the self-energy E" of the internal line (a), and to the tadpole diagrams (b). 



Equations (39) also split to separate equations according to the amount of spin transferee! from the right to the 
left. Let us define 



' A°+ 



+ A_ 1 = 



+A 1 



-A 



A " 

z+ A 

A ( 

L c 

A 



A°° ■ 



Kb 1 

Ao- 1 ' ' 
Ate" 1 , 



(40a) 



(40b) 





FIG. 6: The Feynman diagrams defining the proper anomalous vertex. The hatched triangle represents the proper anomalous 
vertex. 
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for the polar case and 



A ++ \ / A+" 



—a = I A& , "A = A+f I , (41a) 

AT" / V A+" 



ro,i i ! 



-i 



+ A_! =0, +A = 0, (41b) 

q a 1 = a:+, ^A^a-l 1 , 

q A_i = 0, Q A _1 = 0, (41c) 

for the ferromagnetic one. The resulting equation for ™A Q reads as: 

"AJfc) = "A (0)Q + |^ n g {Q) (q, k + q) n g(q, k + q, -q) "A (0)q (42) 

for a = ±1. Where the A, ) vectors are defined similarly as the A vectors in Eqs. (40) and (41). We have used the 
consequence of Eq. (31a), namely that does not depend on its third argument. The solution of Eq. (42) reads as: 



n Aa(k) = n i~\k) »A (0)Q . (43) 

■ — (x 

Similar considerations can be carried on to ™A (fc) with the result: 

n T{k) = n g*~\k) "Af 0) . (44) 



D. Self-energies in the Hartree-Fock approximation 



Until now we have considered the Hartrcc-Fock propagators, which are used in the equation of state. The thermal 
Green's function of the system, defined by Eq. (7) and expressed with the help of the self-energies by the generalized 
Dyson-Beliaev equation (9), are more complicated. To be consistent both with the partial summation made in the 
equation of state and the consequences of symmetry breaking, one has to construct the self-energies from the tadpole 
diagrams in Fig. 5. b) by changing a circle (representing a condensate atom) to an appropriate particle line. The 
self-energies obtained are diagonal in the upper (Roman) indices; moreover this index is in the polar and + in the 
ferromagnetic case, i.e. agrees with that of the condensate. The self-energies and S^T 1 " satisfy by construction 
the relationship: 

-^=^ c o.i = EfitfJ.O) - £-.1(0,0), (45) 

where c = for the polar and c = + for the ferromagnetic case. The Eq. (45) amounts to a generalization of the 
Hugcnholtz-Pines theorem [29] to the spin-1 systems, since 

£§,! = (), so S 1 T(0,0)-S c _ c ia (0,0) = Ai , (46) 

where £' is the self-energy without the subtraction /i. This ensures that the corresponding excitations are gapless as 
required by the general Goldstone theorem due to the breaking of the gauge symmetry. 

As for the other self-energies consistency requires that they should be built by using the same diagrammatic blocks 
and to the same order in them. The resulting Feynman diagrams are depicted in Fig. 7. a) and b). The more detailed 
diagrammatic equation is shown in appendix A. (Note that in the ferromagnetic phase E^j 1 (fc) = except when 
r = s = + due to the structure of the interaction). One expects that S5?(0, 0) = in the ferromagnetic phase and 
£7-T(0, 0) — £+i~i(0 5 0) = in the polar phase. These relationships lead to gapless spin wave spectra as required by 
the Goldstone theorem, since the rotational symmetry in the spin space is also broken, besides the gauge symmetry. 
The spectra to be calculated in the next section will exhibit this property. 
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(b) 



FIG. 7: The Feynman diagrams contributing to the self-energies 



According to Fig. 7. a) the self-energy contributions can be brought into two categories. Those graphs are in 
the first category, which can not be split into two (connected to external lines) by cutting a single interaction line. 
The first four diagrams fall under this category. We call these self-energy contributions proper. The second category 
contains those graphs, which do not exhibit this property. The last diagram in Fig. 7. a) falls under this category. 
We call these latter self-energy contributions improper. Therefore the equation of the normal self-energies read as: 



s;;(fc) = + M[f(fc), with 



E^(fc) = 7T 1 



n + N ei r '(o, k, o) + H r ' s ' (V; s v + v;,';) 



(47a) 
(47b) 
(47c) 



where tilde refers to proper diagrams as in Subsection II C 4. The resulting diagrams for 1 are seen in Fig. 7. b) 
and read as: 



YT hl {k) = 57 s i ; i(fc) + ML\ A {k), with 

£-i,i(*0 = k-'No [e s r ' r s (-k, 0, k) - V r Y s 

ML\ A (k) = H-^-\k)W e d f c (k)A^(k). 



(48a) 
(48b) 
(48c) 



Finally, for future use we remark that the Green's functions can be categorized according to the amount of spin 
they transfer relative to the spin of a condensate atom. For this aim we need the following considerations. Equation 
(9) splits into independent parts according to the amount of spin transfered [19]. For the polar case, let us define the 
following matrices: 



yi.i y l,-l 

Q+- Q++ 

»-i,i 



7(5 



7 5 



°V — 
Zjjg — 



J 7<5 



voo vOO 

^1,1 ^1,-1 

voo' vOO' 

^-1,1 ^-1,-1 

^1,1 ^1,-1 

^-1.1 ^-1,-1 



7(5 



7(5 



(49a) 
(49b) 



with the indices in the upper left corner corresponding to the amount of spin transfered. For the ferromagnetic case 
we define 





r a ++ 

_ C++ 
- y -l,l 


+ G 


= t/i,!, 


Qg 


= Gi.i, 



°V — 
ZjjS — 



75 



V++ V++ 

^1.1 ^1,-1 
^-1.1 ^-1,-1 



7(5 



+ v _ v 00 
— ^1,1, 



(50a) 

(50b) 
(50c) 
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With the use of the above matrices, Eq. (9) splits to independent matrix equations, all of the form 

n g a7 (k) = "e (0) a 7 (fc) + n G ( o )ap (k) n z pa (k) n g^(k), (51) 

with n being 0, + or Q. (Note however, that in the ferromagnetic case for n = + and n = Q anomalous Green's 
functions do not arise, i.e. a = 7 = 1.) 

E. Collective excitations and generalized density correlation functions 

To describe collective excitations one can introduce the correlation functions as follows [19]: 

D nn (k,T) := -(TV[n(k,r)n t (k,0)]), (52a) 
D„(k,r) := -(T T [^(k,r)^(k,0)]), (52b) 
D ++ (k,r) := -(r T [^+(k,r)^(k,0)]), (52c) 
I>2 + (k,T) := -(T T [^(k,r)^ t (k,0)]) ) (52d) 
for k ^ 0, where the density operator: 

n(k) = ^2 a t( k + q) fl r(q), (53a) 
q 

the spin density operator: 

^(k) = Yl a r( k + q)(^)r.aa(q), (53b) 
q 

a non self-adjoint operator generating spin waves: 

^+(k) = ]T 4 (k + q)(F + ) rs a s (q), (53c) 
q 

and another one for quadrupolar spin waves: 

^(k)=^ot(k + q)(^)„o.(q), (53d) 
q 

with F + = F x + iF y . After Fourier transformations these correlation functions can be expressed in the following way 
[19]: 

D nn (k) = ^ 1 j °£(fc) ^ 1 j , (54a) 
D zz (k) = ^ ^ j °^( fc ) ^ °! j ' (54b) 

^++(fc) = 2 ( ! ) + ^( fc ) ( 1 ) < ( 54c ) 

D® + (k) = 4 Q R(k). (54d) 

The matrices on the right hand side satisfy the equations 

n R(k) = h n U{k) + n U{k) n a n £L{k), (55) 
where n C has been defined by Eqs. (35) The polarization function can be written as 

n n(fc) - n n (r) (fc) + n- 1 n g ai {k) n l a {k) o n V(k). (56) 
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Here the first term on the right hand side represents the one-particle line irreducible contributions, which in the 
Hartree-Fock approximation are given by Eq. (33). Furthermore A is defined in the Subsection II C 4 and Q stands 
for the proper Green's function obtained as the solution of Eq. (51) with replacing E by its proper part (given in the 
Hartree-Fock approximation in Eqs. (47b) and (48b)). 
By introducing the so called dielectric function 

"#) = "I - n U(k) n a, (57) 

the frequencies and the damping of collective excitations (after analytic continuation to the second Riemann sheet) 
can be calculated from the equation 

det "|(k,w) = 0, (58) 

for n = {0, +}. For n = Q they are given by the poles of ^11, since ®C_ = 0. 

It is important that the one-particle and collective spectra almost coincide as can be proven by rearrangements of 
the contributions [19]. The spin density relaxation mode in the polar phase, described by the correlation function 
(54b), is an exception, which does not have a counterpart among the one-particle excitations. One can easily convince 
oneself that the general treatment in [19] directly applies to the Hartree-Fock model worked out in the present paper. 



III. APPLICATIONS 



In this section we evaluate the self-consistent Hartree-Fock approximation, discussed in the previous section, to 
parameter values relevant to possible experimental realizations. In most cases we use dimensionless quantities. Let 
us introduce the following dimensionless small parameters: 

e„ := /3 a Nc n , (59a) 
e s :=0oN\c s \, (59b) 
s = \cs\_ = £^ (59c) 

with Pq := l/(fcsTo) the inverse transition temperature and N the total particle number density. The transition 
temperature T is defined as the temperature, where the symmetric phase of the ideal spin-1 Bose gas looses its 
stability, i.e. where 

rtf)F(io) 

with A := A(T = T ) the thermal wavelength (15) of the atom at temperature T . The self-energies are made 
dimensionless with the help of the inverse temperature in the following way: 

:= /3ft££(fc). (61) 

The condensate density and the thermal densities in their dimensionless form read as: 

70 := l3c n N 0l (62a) 



C r ■•= dc n H rr = ^|||^ , (62b) 



where j3 — l/(ksT) the inverse temperature. Automatic summation is not applied to Eq. (62b), the temperature is 
also used in a dimensionless form, with t — T/To and £(s) is the Riemann zeta function. The last equation is obtained 
with the use of Eq. (14). The dependence on the wavenumber is also evaluated in a dimensionless form. For this aim 
the thermal wavelength is used again 

u := |k|A. (63) 

One has to recall at this place that the Hartree-Fock approximation is not applicable in the vicinity of the phase 
transition. Furthermore, it is not valid at very low temperatures where the contributions of Beliaev processes might 
become important like in the case of a scalar condensate [30] . Nonetheless there is a considerable temperature range 
where one can expect that such a Hartree-Fock approach gives quite accurate results. 
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A. Ferromagnetic system 



This system is realized by the gas of 87 Rb atoms, which has been recently trapped and condensed in an all optical 
way [31]. As stated earlier, the condensate spinor is chosen to be £ r = (1,0, 0) T . We consider a unit volume system 
with transition temperature at T = 200 nK (close to that analytically calculated in Ref. [32] for the experimental 
situation realized in Ref. [31]), which by Eq. (60) yields a homogeneous total density N — 10 20 m~ 3 . We take the 
scattering lengths as Oo = 104ae, and a 2 = 102ae ( with clq being the Bohr-radius) according to the calculations 
made in Ref. [33] ( note that their results have some uncertainty, since the exact number of bound states in the triplet 
potential of two 87 Rb atoms is not known exactly). According to Eqs. (59), the values of the dimcnsionless small 
parameters are: e„ = 0.2 and S = 0.0065. 

The length of the condensate is the longest length scale of the system. Taking the length of the condensate as 100 /im 
means a minimal wavenumber fc m in = 10 4 m _1 . For experimental applications local excitations should be considered. 
One can take the length scale of local perturbations to 10 ^m resulting in a typical wavenumber k typ = 10 5 m _1 . 

With the use of Eqs. (13), (59c), (61), (62) and (17) the dimensionless self-energies of the internal lines read as: 



d++ = (1 - 5)7o, 



'ii 

CT 11 



a 



ii 



= -(I - 6)C+ + (1 + S)C - 26C-, 
= -(1 - 3(5)C+ + (1 - 3S)C- + 2(57o. 



(64a) 
(64b) 
(64c) 



Expressing the condensate density as the total density minus the noncondensate densities leads to the equation: 



7o = -j — C + - C - C_ . 



(64d) 



Equations (64) and (62b) form a closed set of nonlinear equations for the densities and self-energies , which can be 
solved numerically. 



4 
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FIG. 8: The characteristic lengths of the ferromagnetic 87 Rb system as functions of the dimensionless temperature t = T/Tq. 



1. Zero spin transfer excitations 



Equations (54a) and (54b) show, that density and spin density collective excitations fall to this category together 
with the one-particle excitations determined by the zero spin transfer Green's function Eq. (50a). According to Eq. 
(31a) the matrix ™0(fci, fe, q) docs not depend on q, moreover it depends only on ki — k\. Using this property and 
Eqs. (47b), (48b), (62) and (17), the dimensionless proper self-energy of spin-transfer can be brought to the form: 

°a Q7 (fc) - /3^V [9++(0, k, 0) - (c + c s )] , (65) 
for all a, 7. The zeroth order contributions of the anomalous vertex vectors (41a) read as: 
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With the help of Eqs. (31a), (31b), (33), (43) and (44) the regular polarization, the anomalous vertex and the 
proper self-energy can be given explicitly. With the use of these and Eqs. (56) and (57) the dielectric function can 
be expressed explicitly. The calculation is straightforward but rather lengthy. According to Eq. (58) one way to 
determine the spectrum of density and spin-density excitations is provided by the equation: 

det °|(k, u) = 0. (67) 

Equation (67) can be directly solved numerically. Analytically solvable approximations can also be obtained by keeping 
the regular terms and the contribution of the nearest pole to the real axis from the Mittag-LefHer representation of the 
spectral function of the bubble graphs occurring [19]. This approximation is valid for a linear excitational spectrum 
when the mean-field characteristic lengths are the largest length scales in the system, i.e. when 3> A, Here 
£n s = ^/-\/4MAo|c„. s | and £' = M/(4-irH 2 N /3). The characteristic lengths of the system are depicted in Fig. 8. It is 
visible that the analytical approximation can be used for t > 0.35. Both the real and imaginary part of the solution 
of Eq. (67) is found to depend linearly on momentum for the relevant range, the solution can be written in the form: 

w = c |k| -ir |k|, (68) 

where c is the speed of the excitation and T is also a positive constant yielding the damping. The numerical and 
analytical solutions are plotted in Fig. 9. These kind of excitations (density and spin density excitations) show similar 
features than the density excitations of the scalar gas [17]. They have linear dispersion curve for small wave number 
with an also linear and small imaginary part. 

c -r 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

* = T/T t = T/T 



FIG. 9: The sound velocity c, and the parameter — T of the 87 Rb ferromagnetic system in units of fesToAo/fi, and as the 
function of t = T/T . 



2. Plus two spin transfer excitations 



As indicated by Eq. (54d) quadrupolar spin waves fall into this category together with the one particle excitations 
given by the +2 spin transfer Green's function Eq. (50c). Spin conservation forbids here the existence of anomalous 
Green's functions. Furthermore since ®C — both the Green's function ®Q and the correlation function ®D are 
proper. Using Eqs. (47b), (62) and (17), the dimensionless proper self-energy of spin-transfer 2 can be brought to the 
form: 

Q a n (k) = /3N Q [ett(0, k, 0) - (c„ + c s )] + (1 - 35)(C_ - C+). (69) 
With the help of Eqs. (28), (31a) and (31b) the 6 matrix for n = Q can be cast to the form: 

9++(0, k, 0) = Cn(l + 6) 

According to Eqs. (51), (61), (69) and (70) the Green's function ®Q is expressed (after the usual analytical continu- 
ation in the frequency) as 

»o% U ) - «s<* U ) = , !(0 _ 1) :\ i , kiU ) : an 
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with u defined by Eq. (63), and il := foui/e^. The poles of the Green's function (71) with the use of Eqs. (69), (70) 
and (62) are determined by the following equation: 



u 2 (n - 1) + (1 - 35)(C+ - C-) - 2<5 7o - (1 + 5)c n Tl(+ + _(u, fi) [u 2 (n - 1) + (1 - 3<5)(C+ - CL) + 7o (l - 5)] = 0. (72) 

Equation (72) can be solved in a numerical way for Q as a function of u. Analytical approximate solution is also 
possible with the help of Eqs. (B3) and (B4). If one is also interested in the damping rates, the small imaginary part 
(B3c) of the bubble graph is also to be considered. The following form is motivated in Appendix B: 

c n 3mU^+ + _(u,n) = e -(^)\. ..), foramft^O, (73) 

with Act = ctJT" — ct^ + . (We have supposed that the imaginary part of SI is very small, therefore 3m U = II".) The 
quantity in the brackets is not important for our considerations here. 




FIG. 10: The gap A and effective mass of the quadrupolar spin wave excitations as functions of the temperature t — T/To. 
The gap is measured in units of fcsTo, while M* is measured in the units of the bare mass of the 87 Rb atom. 

The long wavelength excitations can be written in the form: 

A fik 2 

with A being the energy gap of the quadrupolar excitations, M* their effective mass and 7(|k|) is their wavenumber 
dependent damping rate. The gap A and the effective mass is plotted in Fig. 10. The imaginary part of the 
numerical solution at T = 0.7 To temperature is plotted in Fig. 11. According to Eq. (73) and the numerical 
solution of Eqs. (64), 7(|k|) vanishes exponentially for |k| — > 0. This exponential suppression dominates below the 
characteristic value: 

Uc = l g iT M 073 (75) 

For illustration we have fitted the function 

/(«) = -Cei^T (76) 
to the numerical values in Fig. 11, and found C w 0.089 and u c ss 0.068 in good agreement with (75). 



3. Plus one spin transfer excitations 



As indicated by Eq. (54c) spin waves fall into this category together with the one particle excitations determined 
by the +1 spin transfer Green's function Eq. (50b). Spin conservation forbids here the existence of anomalous Green's 
functions. Using Eqs. (47b), (62) and (17), the dimensionless proper self-energy of spin-transfer 1 can be brought to 
the form: 

+CT n (fc) = I3N Q [e++(0, k, 0) - (c„ + c a )] - (1 - S)C+ + (1 + 5)C Q - 2SC-. (77) 
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FIG. 11: The damping rate of the quadrupolar spin wave excitations at temperature T = 0.7 T . The figure contains the data 
points of the numerical solution and the curve of the fitting function (76). 



The zeroth order contribution for the anomalous vertex vectors Eq. (41b) are 

+ A (0)il = ( J ) , + A (0 ),-i = 0. 



(78) 



Again with the help of Eqs. (31a), (31b), (33), (43) and (44) the regular polarization, the anomalous vertex and the 
proper self-energy can be given explicitly. With the use of these and Eqs. (56) and (57) the dielectric function can 
be expressed. The spectrum of the spin wave excitations is given as the solution of the equation: 

det +|(k, lu) = 0, (79) 
The energy of the excitations is found to be quadratic for small wavenumber (including wavenumber values considered 
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FIG. 12: The effective mass of the spin wave excitations in units of the bare mass of the 87 Rb atom as a function of t = T/Tq. 



in this paper). The spectrum of excitations can be written in the form: 

fik 2 



iv = 



2M+ 



H(|k|), 



(80) 



with M+ the effective mass and 7(|k|) the damping rate. The temperature dependence of the effective mass is plotted 
in Fig. 12, while the damping rate is plotted for T = 0.8 T in Fig. 13. The behavior of the damping rate of the spin 
wave excitations is more complex than that of the quadrupolar spin waves shown in Fig. 11. The reason is that in 
this case two different bubble graphs emerge. One is 14^ with self-energy difference Aa\ — a®® — the other is 

IIq ° with self-energy difference Aa 2 = o^i ~ ^11 ■ ^ s a consequence there are two characteristic wavenumber values, 
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FIG. 13: The damping rate divided by u 2 of the spin wave excitations in units of fcsT/ft at T = 0.8 To. 

where the damping rates suffer a drastic change. At t = 0.8 these dimensionless wavenumbers are: 

Uc.i = 11 — — « 0.0005, (81a) 

u c , 2 = 1(711 1 « 0.05. (81b) 

Since the first wavenumbcr value is much smaller than the minimal accessible wavenumber u Ct i -C Afc m ; n , there is 
some damping of the long wavelength spin wave excitations even, when quadrupolar spin waves are not damped. 

B. Polar system 

This system is realized by the ultracold gas of 23 Na atoms, when the magnetic field is sufficiently small [3, 4]. 
We consider a pure condensate and take its spinor as ( r = (0, 1,0) T . The system considered is also of unit volume 
with periodic boundary conditions and a transition temperature of T = 2 fiK. The total particle density calculated 
with Eq. (60) is N = 4.6 • 10 20 m -3 . We take the scattering lengths of the sodium atoms from Ref. [34] to be 
a = 50a# and a 2 = 55as. According to Eqs. (59), the values of the dimensionless small parameters are: e„ = 0.17 
and S = 0.031. To determine the interesting wavenumber range, we also take the condensate size to 100 /xm, from 
which the minimal wavenumber is fc m ; n = 10 4 m _1 . The length scale of local perturbations is chosen to correspond to 
a typical wavenumber value fc typ = 10 5 m _1 . 

With the use of Eqs. (13), (61), (62) and Eq. (20), and the (+, — ) symmetry of the polar case [namely that 
(T]Y~ = the dimensionless self-energies of the internal lines read as: 

3% - 70, (82a) 
5=++ - (l-S)(C+ -Co) + 5 lQ . (82b) 

Expressing again the condensate density as the total density minus the noncondensate densities and using the sym- 
metry mentioned leads to the final equation: 

7o = y - 2C+ - C . (82c) 

Equations (82) and (62b) are closed and numerically solvable for the densities and self-energies in the polar phase. 

The characteristic lengths of the polar system are plotted in Fig. 14. As can be seen from the figure, the approxi- 
mation is usable roughly for t > 0.3 (when the largest characteristic lengths are the mean field correlation lengths). 



Zero spin transfer excitations 



As in the case of the ferromagnetic system, these type of excitations are the particle density and spin density 
excitations defined as the poles of the correlation functions (52a) and (52b) [See Eqs. (54a) and (54b)]. Since the 
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FIG. 14: The characteristic lengths of the polar 23 Na system as a function of temperature t = T/Tq. 



matrix does not depend on the momentum transferred [See the remark after Eqs. (31)], the proper self-energy of 
spin-transfer for the polar system can be cast to the form (by using Eqs. (47b), (48b), (62) and (20)): 

oF a7 (fc)= J 9JV [e§8(0,fc,0)-c n ], (83) 



for all a, 7. The zeroth order contributions of the anomalous vertex vectors (40a) read as: 

°A(0), Q - V^Vo (Vj ■ (84) 

With a lengthy but trivial calculation, using Eqs. (31a), (31b), (33), (43), (44), (56) and (57), one can give the explicitc 
form of the dielectric function °e(k, iuj n ). Analytically continuing it in the usual way in the u> variable one determines 
the elementary excitations by 

dct °£(k,w) = 0. (85) 



It can be seen from the detailed calculation, that the left hand side of Eq. (85) for the polar case is the product 
of two terms. One term determines the particle density excitations and the other term determines the spin density 
excitations. This behavior is not limited to this approximation. As shown generally in Ref. [19], the particle density 
and spin density correlation functions D nn and D zz , respectively, have different excitation spectra in the polar case, 
even if they both belong to spin transfer zero. 

The spin density excitations are the solutions of the following equation: 

l-(l + 3<5)c„n++ ++ (n,O) = 0, (86) 

with O = fiw/ek and u = |k|A. This is an overdamped mode, its frequency is purely imaginary with linear wavenumber 
dependence: 

w = -ir a |k|, (87) 

where r s is a positive parameter with the dimensionality of velocity. The numerical and analytical solution is plotted 
in Fig. 15. The analytical solution is obtained from the Mittag-LefHer expansion [16, 19] of the spectral function of 
the bubble graph. 

The particle density excitations show the same character as their ferromagnetic counterparts (or as the density 
excitations in a scalar gas). Namely they have a linear energy spectrum and linear damping rate in the small 
wavenumber limit. The numerical solution and the analytical approximation from the Mittag-LefHer series [16, 19] is 
plotted in Fig. 16, where c and T are defined as in Eq. (68). 
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t = T/T 

FIG. 15: The relaxation rate parameter — F s of the spin density wave of the 23 Na (polar) system in units of kBToXo/ft, and as 
the function of t = T/T . 
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FIG. 16: The sound velocity c, and the parameter — T of the particle density excitations of the Na (polar) system in units of 
fesToAo/fi, and as the function of t — T/To. 



2. Plus one spin transfer excitations 



The +1 spin transfer excitations are spin waves, as shown by Eq. (54c). Using Eqs. (47b), (62) and (20), the 
dimcnsionless proper self-energy of spin-transfer 1 can be brought to the form: 

+a 11 (k) = (3N [e°°_(0,k,0) - c n ] + (1 - 5)(C+ - C ), (88) 

+5 1) _ 1 (fc)=/37V [et°(0,fc,-/e)- Cs ], (89) 

+ oF_i,i(fc) = /3N [e° Q ±(-k,0,k) - c s ], (90) 

= /3iVo[e+ + (-fc,0,0) - c n ] + (1 - S)(C + - Co). (91) 
The zeroth order contribution for the anomalous vertex vectors Eq. (40b) are 



+ A (0)il = ( J ) , + A (0); _ 1 = v^Vo ( I ) • 



(92) 



Again with the help of Eqs. (31a), (31b), (33), (43) and (44) the regular polarization, the anomalous vertex and the 
proper self-energy can be given explicitly. With the use of these and Eqs. (56) and (57) the dielectric function can 
be expressed. The spin wave excitations of the polar state are given by the equation: 

det +|(k, lu) = 0. (93) 

Since the n°^ +0 and n^ _ bubbles appear in Eq. (93), the analytical solution is also obtained from the asymptotic 
series expansion (B4) of the spectral function of the bubble graph. 
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In the polar case t?^ = <7-jV~, therefore the self energy difference of the internal lines for the two bubbles are equal, 
and as a consequence the damping rate of the spin wave excitations will be exponentially suppressed but with one 
characteristic wavenumber value: u c = — o : J t 1 + |/2. The spectrum of spin wave excitations in the polar case can 
be written in the form: 

c 

0.12 
0.1 
0.08 
0.06 
0.04 
0.02 
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t = T/T 

FIG. 17: The velocity of the spin-wave excitations of the 23 Na (polar) system in units of kBTo\o/ft, and as the function of 
t = T/T . 

u) = c|k| +n(|k|). (94) 

The numerical and analytical solutions for the sound velocity of the spin wave excitations are to be found in Fig. 
17. The damping rate is plotted for T = 0.8 T in Fig. 18. 
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FIG. 18: The damping rate of the spin wave excitations at temperature T = 0.8 T for the 23 Na (polar) system in units of 
(k\) 2 /(h[3), and as the function of u — |k|A. 

As a summary we compare the frequencies and damping rates of the various excitations of both the ferromagnetic 
and polar cases for a given momentum value, namely fc = fc typ = 10 5 m _1 . In Fig. 19. the excitation frequencies 
and damping rates of the collective excitations of the ferromagnetic system are plotted. It can be seen, that for 
this wavenumber value, density excitations (density fluctuations) cost much more energy and are relatively short- 
lived compared to spin waves and quadrupolar spin waves. These latter excitations correspond to transverse spin 
fluctuations of the system with constant density. In Fig. 20. the frequencies and damping rates of the excitations of 
the polar system are shown. Density fluctuations cost more energy than transverse spin fluctuations (in the polar case, 
we mean transverse as orthogonal to the spin quantization axis) in this case too and are also shorter lived than the 
latter ones. The shortest living excitations in the polar case however are the non-propagating spin density excitations, 
which belong to longitudinal spin fluctuations. 
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FIG. 19: The frequencies (e = ^Ktu)) and damping rates (7 = 3mw) as functions of dimensionless temperature t = T/To of 
the various collective excitations for k = k typ = 10 5 m _1 of the ferromagnetic system. Both quantities are measured in units of 
k B T Q /h. 
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FIG. 20: The frequencies (e = Dlcco) and damping rates (7 = 3mw) as functions of dimensionless temperature t = T/To of the 
various collective excitations for k = k typ = 10 5 m _1 of the polar system. Both quantities are measured in units of feT /fi. 



IV. DISCUSSION 



In this paper we have generalized the self-consistent Hartree-Fock theory to the symmetry breaking phase of the 
spin-1 Bose gas. We have started from the Hartree-Fock equation of state of the ferromagnetic and polar condensed 
phases. The Green's functions and correlation functions have been given, in accordance with the self-consistent 
equation of state and with the dielectric formalism, from which the one-particle and collective excitations of the 
system have been determined. The results have been applied to the cases of 23 Na (polar) and 87 Rb (ferromagnetic) 
gases. 

Both in the ferromagnetic and polar cases the behavior of particle density excitations, belonging to spin transfer 0, 
has been found similar to that of the particle density waves of scalar Bose gases, where these are the solely existing 
collective modes. In particular both their frequency and damping rate depend linearly on wavenumber, see Eq. (68), 
Figs. 9. and 16., and Ref. [12, 16, 17, 19] in the intermediate temperature region. These collective excitations 
agree with one-particle excitations belonging to spin transfer due to hybridization caused by the condensate. In 
the polar case also for spin transfer another collective mode (independent of particle density waves) exists, which is 
overdamped (Fig. 15.) and represents spin relaxation. This collective mode is related to the pole of the correlation 
function (52b). Furthermore, it is absent from the one-particle excitations [19], which is the consequence of the + <-> — 
spin inversion symmetry of the polar phase and is independent of the approximation. Spin waves are described by the 
correlation functions (52c) and are collective excitations with spin transfer 1. They hybridize with the one-particle 
excitations with the same spin transfer in both the ferromagnetic and polar cases. It has been found that with taking 
into account the exchange processes the properties of these excitations calculated in Ref. [19] essentially change. In 
the ferromagnetic case, though the excitation energies continue to show a quadratically starting energy spectrum, 
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the behavior of the effective mass as a function of temperature (Fig. 12.) differs qualitatively from that of the same 
excitations calculated in the Hartree approximation, which contains only direct interaction processes. Note, that 
in the Bogoliubov approximation the effective mass has been the bare one [8, 9]. In the polar case the excitation 
frequencies depend linearly on wavenumber with a velocity depicted on Fig. 17. This feature has been found already 
in the Bogoliubov and the Hartree approximations [8, 9, 19]. In these cases, however, the interaction parameter c„ 
does not show up in the velocity of the spin wave. It is no more the case in the present Hartree-Fock approach and 
the solution drawn on Fig. 17. includes it. The imaginary part of the spectrum, i.e. the damping rate is, however, 
exponentially suppressed (Fig. 18.), differing qualitatively from the damping rate in the Hartree approximation, 
which is linearly dependent on wavenumber in the relevant temperature region. The reason of the observed significant 
difference is that, while in the Hartree approximation all three self-energies belonging to internal lines with different 
spin projection are zero, in the Hartree-Fock approximation only the self-energies of the and propagators 
are equal in the polar case. This latter coincidence follows from the + <-» — symmetry of the polar state. The lift 
of the degeneration results in the appearance of the characteristic wavenumber under which the damping rate of the 
excitations are suppressed. (In the ferromagnetic case all three internal lines have different self-energies, therefore 
the damping rate of both the spin wave and quadrupolar spin wave are suppressed for small momentum values, see 
Figs. 11. and 13.) This can be regarded as a threshold effect. One expects that in the intermediate temperature 
region its occurrence is not restricted to the Hartree-Fock approximation, but it is generally present though can be 
less pronounced. The reason is that in this temperature region scattering processes at wavenumber greater than 
the inverse of the Bogoliubov coherence length make the important contributions and the use of the Hartree-Fock 
propagator for the internal lines is justified. Collective excitations with spin transfer 2 (quadrupolar excitations) are 
propagating only in the ferromagnetic case. The frequency of these excitations start with a gap (since these are non 
Goldstone excitations) and is parabolic (see Fig. 10.). Its effective mass is close to that of the ferromagnetic spin 
waves, the difference is proportional to the smaller coupling constant c s . The existence of a gap and the quadratic 
dispersion has been found already in the Bogoliubov and Hartree approximations [8, 9, 19]. 

Some remarks concerning hybridization of one-particle and collective excitations are in order here. We have used 
this terminology in the sense as introduced in case of a scalar condensate. Namely, that due to the presence of the 
condensate one can see intermediate states corresponding to collective excitations in the perturbation series of the 
one-particle propagators and similarly, one can identify one-particle intermediate states in the perturbation scries 
of the correlation functions describing collective modes. As a result the one-particle and the collective excitations 
spectra coincide after the necessary rearrangement of the perturbation series, carried out within the framework of 
the so called dielectric formalism (see for a review Ref. [13]). In case of the systems with a spinor condensate there 
are exceptions when such hybridization docs not take place [19]. One example is the spin density relaxation in the 
polar phase described by the correlation function D zz (52b), whose perturbation series does not contain one-particle- 
like intermediate states (see Eqs. (86), (87)). Another example is provided by the quadrupolar excitations in the 
ferromagnetic phase. Though in the perturbation series D® + (52d) one-particle like intermediate states occur, it is 
not true vice versa since the Green's function is proper in this case (see III A 2. chapter of the Applications). As a 
result the spectra of the correlation functions (52) is more complete than those of the Green's functions. A further 
advantage of using the correlation functions is that the physical nature of the excitation calculated can be seen from 
the definition of the relevant correlation function. By these reasons we also have put emphasis on them in the Section 
devoted to the applications. 

We should also comment on the region of validity of the self-consistent Hartree-Fock approximation. This approx- 
imation can not be used very close to zero temperature. In this regime other processes give significant contributions 
to the self-energies (Beliaev processes). The Hartree-Fock approximation is expected to provide quantitatively good 
results when the mean-field coherence lengths are the longest length scales in the system. It can be seen in Figs. 
8. and 14. that this is the case above T w 0.3Xo- In this model approximation it is found that the transition is a 
first order one for both the ferromagnetic and polar cases, such as for the scalar gas [17]. On the other hand the 
transition is proved to be second order for the polar case and weakly first order for the ferromagnetic case in the 
self-consistent Hartree approximation [19]. Since both approximations represent improved mean-field approaches they 
are not applicable near the transition point, where a more sophisticated treatment is needed, which can also decide 
upon the order of transition. 

V. ACKNOWLEDGEMENT 

The present work has been partially supported by the Hungarian Research National Foundation under grant No. 
OTKA T029552. 



24 



APPENDIX A: THE SELF-ENERGIES IN THE HARTREE-FOCK APPROXIMATION 

In this appendix we show how the self-energies X^f and T 1 c 2 ll (with c = for the polar and c = + for the 
ferromagnetic phase), depicted in general in Fig. 7., can be obtained by changing all condensate circles to external 
lines, with the same spin projection, in the tadpole diagrams Fig. 5. b). The result of the above process is illustrated 
in Fig. 21. With the help of the diagramatic equation in Fig. 5. a) one can replace the self-energy contribution of 
the internal line by the simpler form depicted in Fig. 1. Also using the anomalous vertex vectors, plotted in Fig. 6 
one arrives at the self-energies as illustrated in Fig. 7. 
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FIG. 21: The structure of the self-energies Efj and Effi i 



APPENDIX B: THE ANALYTICAL APPROXIMATION OF THE BUBBLE GRAPHS WITH 

DIFFERENT SELF-ENERGIES 

In this appendix we outline the approximation of the bubble graphs with different self-energies of the propagators. 
The contribution of the analytically continued bubble graph, with internal propagators carrying spin s and spin r, 
forward and backward, respectively, can be cast to the form [19]: 



C» n (0)r«(«> n ) = 



12r(3/2)C(3/2)w 



rVo rVi 

/ Rarr( z )dz- / Rzss(z)d 

J oo J oo 



(Bla) 



25 



with the notations 



y = l-[u 2 (Cl-l)-Aa}, (Bib) 

?/i = ^-[u 2 (n + 1) - Act], and (Blc) 
Act = ct^-ct[[. (Bid) 



The function 



f°° 2x 1 

RAZ) = L x e* 2 +*-lx-z dX > (f ° r 3mz> °) ( B2 ) 

(analytical in the upper half-plane) is to be continued to the whole complex plane. We are interested in the 3m ft < 
values, since the elementary excitations have finite lifetimes. Therefore z also has a negative imaginary part, therefore 
the function reads as: 

J? s (z) = R~ (z) + &i (z) , with (B3a) 
RUZ) = loo e^-lx-z dX (B3b) 

4? 7T 7 

R'i{z) = ez2+ ~ a _ r (for 3mz<0). (B3c) 

In the long wavelength limit, i.e. when 2m <C |Act|, the upper limits of the integrals in Eq. (Bla) diverge, therefore 
an asymptotic series expansion can be used to represent the integral Eq. (B3b): 

jv m 2 ~ r(n + 3/2)F(n + 3/2,CT) 

n=0 

Let us define the quantities n^ rg and n^. g , with Eq. (Bl) by replacing R with R' , and R", respectively. One 
obtains 



2irie n Vt 1 - e^o +<T i 

12r(3/2)C(3/2)u n l - e-^ii ' 



C „n' (0 ^( U , O) = 1oWo/o ;,; o/ ^, In - _ 24 _ 3as ■ (B5) 



In the same limit, i.e. when 2u <C |Act|, one can further approximate c n H^ rs with a result: 



27rie„v / i 
12r(3/2)C(3/2)iT 



w 2 (fi - l) 2 - 2AS(fl - 1) - 4ctJ| 



exp 



4 

w 2 (f2 + l) 2 - 2Aa{Q + 1) - 4ct[[ 



• (B6) 



Taking the real part of f2 positive, and with the assumption that its imaginary part is small one can see that the 
quantity in the square brackets is not important for the considerations made in Sec. III. 

We have considered the analytical approximation of the retarded part of the bubble graph's contribution with dif- 
fering self-energies in the forward and backward propagating lines, i.e. with nonzero spin transfer. The approximation 
for such bubble graphs with the same self-energies in their internal lines, i.e. those occurring in zero spin transfer 
processes, are obtained by a Mittag-Leffler expansion, discussed in detail e.g. in Ref. [16]. 
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